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Abstract. 

The weak-coupling expansion of the QCD free energy is known to order g^, logg s , however, the resulting series is poorly 
convergent at phenomenologically relevant temperatures. In this proceedings, I discuss hard-thermal-loop perturbation theory 
(HTLpt) which is a gauge-invariant reorganization of the perturbative expansion for gauge theories. I review a recent NNLO 
HTLpt calculation of QCD thermodynamic functions. I show that the NNLO HTLpt results are consistent with lattice data 
down to temperatures ~ 2T C . 
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INTRODUCTION 

The current generation of ultrarelativistic heavy-ion col- 
lision experiments should exceed the energy density nec- 
essary for the formation of a quark-gluon plasma. Initial 
temperatures of RHIC are up to twice the QCD critical 
temperature, T c ^170 MeV. The strong coupling con- 
stant at these initial temperatures is approximately g s ~ 2 
or a s = g^/4n ~ 0.3, which is some intermediate value, 
neither infinitesimally small nor infinitely large. Theo- 
retically, one expected that this state of matter could be 
described in terms of weakly interacting quasiparticles; 
however, data from RHIC suggested that the state of mat- 
ter created there behaved more like a strongly coupled 
fluid with a small viscosity [1]. This has inspired work 
on strongly-coupled formalisms. However, some observ- 
ables such as jet quenching [2] and elliptic flow [3] can 
also be described using perturbative methods and so it 
is difficult to judge whether the plasma is strongly or 
weakly coupled based only on RHIC data. The initial 
temperatures of the upcoming experiments at LHC are 
expected up to 4 — 6 T c and due to asymptotic freedom of 
QCD, this corresponds to a smaller coupling constant. A 
key question is then whether the matter generated can be 
described in terms of weakly interacting quasiparticles at 
these higher temperatures. 

The weak-coupling expansion of the QCD free energy 
is known up to order g^\ogg s [6, 7]. Unfortunately, the 
resulting series shows no sign of convergence at phe- 
nomenologically relevant temperatures. There are sev- 
eral ways of reorganizing the perturbative series at finite 
temperature [8] and they are all based on a quasiparti- 
cle picture where one is perturbing about an ideal gas of 
massive quasiparticles, rather than that of massless par- 
ticles. In the following I will discuss recent advances in 
the application of hard-thermal-loop perturbation theory 
(HTLpt). 



HARD-THERMAL-LOOP 
PERTURBATION THEORY 

Hard-thermal-loop perturbation theory is a gauge- 
invariant extension of screened perturbation theory [9]. 
The basic idea of the technique is to add and subtract 
an effective mass term from the bare Lagrangian, and 
to associate the added piece with the free Lagrangian 
and the subtracted piece with the interactions. In gauge 
theories, however, simply adding and subtracting a local 
mass term violates gauge invariance [10]. Instead one 
adds and subtracts an HTL improvement term, which 
dresses the propagators and vertices self-consistently 
so that the reorganization is manifestly gauge invari- 
ant [11]. HTLpt has recently been pushed to NNLO 
and the details of the formalism and calculations are 
presented in Refs. [12, 13, 14]. Here only a few selected 
results from QCD [14] are reviewed. 

With rescaled dimensionless parameters m D j q = 
m D i q / (2nT) and fx = h/(2jcT), the renormalized 
NNLO thermodynamic potential for QCD with Nf 
flavors and N c colors reads 
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FIGURE 1. Comparison of LO, NLO, and NNLO predictions for the scaled pressure for Nf - 
panel) with lattice data from Bazavov et al. [4] and Borsanyi et al. [5]. See main text for details. 
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In order to complete a calculation, a prescription is 
required to determine the mass parameters mg and m q . 
Here the Debye mass is set to the mass parameter of 
three-dimensional electric QCD (EQCD) [7], i.e. mo — 
m.E- In Ref. [7], it was calculated to NLO giving 
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The quark mass is set to w ? = 0. 

In Fig. 1, 1 show the scaled QCD pressure for Nf = 3 
(left panel) and Nf — 4 (right panel) as a function of T. 
The results at LO, NLO, and NNLO use the BN mass 
given by Eq. (2) as well as m q — 0. For the strong cou- 
pling constant a s , three-loop running [15] with = 



344 MeV [16] is used here. The bands correspond to 
varying the renormalization scale ji by a factor of 2 
around ji = 2kT which are the central lines. 

The lattice data from the Wuppertal-Budapest collab- 
oration uses the stout action. Since their results show 
essentially no dependence on the lattice spacings (it is 
smaller than the statistical errors), they provide a contin- 
uum estimate by averaging the trace anomaly measured 
using their two smallest lattice spacings corresponding 
to N t = 8 and N T = 10 [5]. Using standard lattice tech- 
niques, the continuum-estimated pressure is computed 
from an integral of the trace anomaly. The lattice data 
from the hotQCD collaboration are their N T = 8 results 
using both the asqtad and p4 actions [4]. The hotQCD 
results have not been continuum extrapolated and the er- 
ror bars correspond to only statistical errors and do not 
factor in the systematic error associated with the calcula- 
tion which, for the pressure, is estimated by the hotQCD 
collaboration to be between 5 - 10%. 

As can be seen from Fig. 1 the successive HTLpt 
approximations represent an improvement over that of a 
naive weak-coupling expansion; however, as in the pure- 
glue case [13], the NNLO result represents a significant 
correction to the LO and NLO results. That being said the 
NNLO HTLpt result agrees quite well with the available 
lattice data down to temperatures on the order of 2 T c ~ 
340 MeV for both N f = 3 and Nf = 4. 

In Fig. 2, I show the NNLO approximation to the 
scaled QCD trace anomaly as a function of T for Nf = 3 
(left panel) and Nf = 4 (right panel). The left panel 
shows data from both the Wuppertal-Budapest collabora- 
tion and the hotQCD collaboration taken from the same 
data sets displayed in Fig. 1. In the case of the hotQCD, 
the results for the trace anomaly using the p4 action show 
large lattice size affects at all temperatures shown and the 
asqtad results for the trace anomaly show large lattice 
size effects for T ~ 200 MeV. The right panel displays a 




parameterization (solid blue curve) of the trace anomaly 
for Nf = 4 published by the Wuppertal-Budapest collab- 
oration [5] since the individual data points were not pub- 
lished. Both panels show very good agreement with the 
available lattice data down to temperatures on the order 
of T ~ 2 T c , Note that due to the massless quark descrip- 
tion in HTLpt, deviations are expected from the lattice 
data for T< 414 Me V. 



CONCLUSIONS AND OUTLOOK 

In this proceedings, I briefly reviewed recent NNLO 
results for the QCD thermodynamics using HTLpt. From 
comparison with lattice data for Nf E {3,4}, it has been 
shown that HTLpt is consistent with available lattice data 
down to T ~ 2 T c for the pressure and the trace anomaly. 

In closing, I emphasize that HTLpt provides a gauge 
invariant reorganization of perturbation theory for cal- 
culating static and dynamic quantities in thermal field 
theory. Given the good agreement with lattice data for 
thermodynamics, it would be interesting to apply HTLpt 
to the calculation of real-time quantities at temperatures 
that are relevant for LHC. 
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